Friction at Atomic-Scale Surface Steps: Experiment and Theory 
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Experiments performed by friction force microscopy at atomic-scale surface steps on graphite, 
M0S2, and NaCl in ambient conditions are presented. Both step-down and step- up scans exhibit 
higher frictional forces at the edge, but distinguish in their load dependence: While the additional 
frictional force due to the step edge increases linearly with load if the tip has to jump a step up, 
it remains constant for downward jumps. This phenomena represents a universal effect that can be 
explained in terms of a modified Prandtl-Tomlinson model featuring a Schoebel-Ehrlich barrier at 
steps. 

PACS numbers: 68.35.Af 62.20.Qp 07.79.Sp 68.37.Ps 



Tribology - the science of friction, wear, and lubrica- 
tion - has impact on many fields of science and tech- 
nology. Consequently, it has been the subject of intense 
research during the last centuries [J. With the advent 
of new experimental techniques such as friction force mi- 
croscopy (FFM) ,2], the study of frictional phenomena 
at the atomic scale became accessible to tribologists, and 
the field of nanotribology has been established since (see, 
e.g., Refs.ifl g H). 

The basic paradigm of nanotribological research is that 
the frictional behavior of a single asperity contact needs 
to be clarified in order to better understand friction in 
complex macroscopic systems. Naturally, this makes the 
friction force microscope a tool of choice for nanotribol- 
ogy. So far, most experimental FFM studies designed to 
elucidate the atomic-scale principles of friction focused 
on the frictional behavior of atomically flat surfaces (see, 
e.g. , Refs.0, S, H EH EH, El E3) . The stick-slip phenom- 
ena observed in these experiments can be understood in 
the framework of the well established and surprisingly 
simple Prandtl-Tomlinson (PT) model (13 . 15 1, which 



is sometimes also referred to as independent oscillator 
model 3, 17 1 . An extension of this model, obtained by 
including thermal activation processes, helped to under- 
stand the velocity and temperature dependence of fric- 
tion [II, El III . 

Although these studies enabled valuable insight into 
the origin of atomic-scale friction, atomically flat sur- 
faces represent a simplified model case, since any truly 
advanced model must include roughness. From this 
point of view, it is surprising that only very few studies 
[21 , 2^, 23|, 24 , 25| focusing on the analysis of friction at 
atomic-scale surface steps exist [26| . All of them reported 
increased frictional forces at step edges compared to the 
value found on atomically flat terraces, which has poten- 
tially far-reaching implications on the friction observed 
in macroscopic systems. Meyer et al. [24[ attributed 
this effect to the influence of the Schwoebel-Ehrlich bar- 
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FIG. 1: (Color online) The frictional behavior of graphite at 
atomic-scale surface steps, (a) Top panel: Color-coded top 
view of the topography; image size: 1.6 /im x 0.4 /im. Middle 
panel: Topographical line section along the arrow indicated 
in the top panel. Individual graphene layers are represented 
by grey rectangles for illustration purposes, revealing that a 
double, a fivefold, and a single step are encountered. Bot- 
tom panel: Frictional forces recorded during a left-right scan 
(red) and the corresponding right-left scan (blue) along the 
same line with an external load of Ji oa d = 6.8 nN. An higher 
frictional increase of the frictional forces for upward than for 
downward scans is obvious (black arrows), (b) Plots of the 
frictional increase (defined as the difference between the max- 
imum frictional force at the step edge and the friction en- 
countered on the terrace) observed at the three different step 
edges as a function of the load. For upward scans (top panel), 
the frictional increase grows linearly with load, while it is 
constant for downward scans (bottom panel). Parameters: 
c z = 0.073 N/m, c x = 20.6 N/m, v M = 6/im/s. 



rier [27J, [28| present at atomic step edges. Analyzing 



the load dependence of an SiaN4 tip on graphite in ul- 
trahigh vacuum, Muller et al. |25[ found a directional 
dependence: Despite being much larger compared to the 
on-terrace value, frictional forces were load-independent 
when jumping a step down ("downwards scan"). In con- 
trast, they increased linearly with load if the tip had to 
move a step up ("upwards scan"). 

The results presented in this letter demonstrate that 
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FIG. 2: (Color online) The frictional increase observed at a 
double step on M0S2 (a) and a single step on NaCl (b). On 
both materials, the additional frictional forces caused by the 
step edges increase linearly with load for upward scans while 
they are independent of the actual loading force for down- 
ward scans. Note that in agreement with a previous report 
by Meyer and Amer [2ll |. we observe no change of friction 
at the edge compared to the on-terrace value for downward 
scans on NaCl, which leads to a vanishing value for the fric- 
tional increase. Parameters: c z = 0.065 N/m, c x — 18.2 N/m, 
v M = 2^m/s (a), and c z = 0.063 N/m, c x = 17.8 N/m, 
i>m = 2 (jm/s (b). 



the findings by Miiller et al. are not restricted to graphite 
surfaces under vacuum environment, but appear to be 
of general significance. Our investigations were per- 
formed on freshly cleaved graphite(OOOl), MoS2(001), 
and NaCl(OOl) surfaces using rectangular silicon can- 
tilevers (ContGD, BudgetSensors) in a commercial fric- 
tion force microscope (MultiMode AFM with Nanoscope 
Ilia electronics by Veeco Instruments, Inc.) operated un- 
der ambient conditions. The normal and lateral spring 
constants c z and c x as well as friction and load were de- 
termined applying the calibration procedures described 
in Ref.|U 

Figure [T] summarizes the data obtained on graphite, 
where we measured the frictional forces at three differ- 
ent step edges of one, two, and five graphene layers height 
within the same scan. Maximum frictional forces encoun- 
tered at a step edge were always higher for upwards scans 
than for downward scans [cf. Fig.QTa), bottom panel, for 
illustration] . In order to adequately study the additional 
contribution of the step edge to the overall friction sys- 
tematically as a function of the externally applied load 
-Fload, we always plot in the following the difference be- 
tween the frictional force needed to overcome the surface 
step and the frictional forces on the terraces, which we 
refer to as "frictional increase" . Our analysis of this ad- 
ditional, step edge-induced component revealed that it 
increases linearly for upward scans, while it remains con- 
stant for downward scans [Fig.QJb)]. Note that we re- 
cover the same qualitative relationships for all three step 
heights, even though higher step heights show higher ab- 
solute increases. Control experiments carried out on a 
double step on M0S2 [Fig.^a)] and a single step-edge on 
NaCl(OOl) [Fig. 0(b)] reveal the same behavior. 




FIG. 3: A schematic of the modified Prandtl-Tomlinson 
model describing the friction at atomic-scale surface steps 
(not to scale). x% represents the position of the tip, which 
is connected via a spring with spring constant c x to the body 
M. For sliding, the body M is moved along the x direction 
while the tip interacts with the tip-sample potential V(xt) 
(thick solid line). If x t = im, the spring is in its equilibrium 
position. The tip movement in the interaction potential at 
the step edge is indicated by arrows. For details, see text. 

Combining the previously published reports with our 
experimental findings suggests that the described load 
dependence of friction at step edges is likely to be a 
general phenomena and independent of both the sam- 
ple material as well as the specifics of the environment 
(vacuum or air). However, despite this fairly far-reaching 
suspected range of validity and its consequential impact 
on macroscopic friction, this effect has not been the sub- 
ject of thorough theoretical analysis yet. To fill this gap, 
we propose an extended Prandtl-Tomlinson model that 
includes an explicit description of the tip-sample inter- 
action at atomic-scale surface steps and, as a result, cor- 
rectly reproduces the load dependence at atomic-scale 
surface steps. 

A schematic representation of the model is shown in 
Fig-El A point-like tip is elastically coupled to a main 
body M with a spring possessing a spring constant c x in 
x direction and interacts with the sample surface via an 
interaction potential Vt s {xt, z^), where xt and z t reflect 
the actual position of the tip. The body M experiences a 
constant loading force .Fioad while it is scanned over the 
sample surface with a velocity vm- The path of the tip 
can be calculated from the equations of motion [30( 

m x x t = c x (v M t - x t ) - -J x Xt (la) 
dx t 

ip dVts ■ /,u 

m z z t = -i^oad ^ IzZt , (lb) 

OZt 

where m x , m 2 are the effective masses of the system and 
7x> 7z damping constants. 

Within the PT model, the point-like tip represents the 
average of the actual tip-sample contact, which might 
include several dozens or even hundreds of atoms. The 
effect of the finite contact size is then mapped into the 
specific choice of the interaction potential for calculation 
convenience. In the classical PT model, this interaction 
is typically approximated by a simple sinusoidal term. To 
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include the effect of scanning over a step edge, we extend 
this approach by explicitly introducing the tip-sample in- 
teraction at surface steps into V tB (x t , z t ). Due to the lack 
of a suitable analytical description, we used a numerical 
approach as described below, which ultimately allowed to 
recover all experimentally observed characteristics from 
the simulations. 

The two-dimensional tip-sample interaction 
potential Vt s (xt,z t ) is computed by the sum- 
mation of individual Lenard-Jones potentials 
v ts = Eili E o {(ro/n) 12 - 2(r /r;) 6 ) where n rep- 
resents the distance between the point-like tip and the 
ith surface atom. The parameters E and r$ describe 
the binding energy and the equilibrium distance, re- 
spectively. Figure Ufa) shows a color-coded density 
plot of a tip-sample interaction potential calculated 
for a hexagonal structure with an atomic distance of 
a = 0.3 nm. By introducing this potential into the 
equation of motion |T]) , we can compute the path of the 
tip on the sample surface, obtaining the lateral force 
from F x = c x (x M -Xt). 

Lateral force curves calculated in this way are plotted 
in Fig.[4jb) as a function of the loading force. Left and 
right from the step edge, F x (xm) exhibits for fixed loads 
the typical saw tooth-like shape with the periodicity of 
the atomic lattice expected from the classical PT model. 
This behavior occurs because the condition 
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is fulfilled Q and reproduces the tip behavior observed in 
actual FFM experiments on atomically flat terraces @, 0, 
H, B Ell EL E2] ■ The significantly altered tip movement 
in the immediate vicinity of the surface step, however, 
requires a more careful consideration. 

For analysis, we start by computing the actual z- 
position of the tip from the stability condition i<i oa d = 
— 9Vt s /(9z t . Calculating the tip-sample potential at these 
positions, we obtain the potential the tip experiences 
for different loading forces [Fig.QJc)]. Caused by a re- 
duction of the atomic coordination at the step edge, a 
step - induced potential barrier (Schwoebel-Ehrlich barrier 
271 . |28| ) is obtained. Conversely, increasing coordination 
leads to a potential minimum at the bottom of the step. 
Left and right from the surface step, the potential shows 
the sinusoidal shape assumed in the classical PT model. 

The lateral force curves plotted in Fig.UJb) can now 
be understood by imagining that the tip moves in this 
potential landscape during scanning. For a downward 
scan [left to right in Fig.QJc)], the tip needs three jumps 
to overcome the surface step [dashed arrows in Figs. [3] 
and|4jc)]. First, it jumps over the Schwoebel-Ehrlich 
barrier into the minimum at the bottom of the step edge. 
The mechanism behind this movement is the same as for 
the earlier discussed stick-slip movement on the terraces: 
The tip sticks at the left slope of the maximum until 
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FIG. 4: (Color online) (a) Density plot of the tip-sample in- 
teraction potential assumed in the simulations, featuring a 
barrier at the top of the step edge and a distinguished poten- 
tial minimum at its bottom (dark blue), (b) The frictional 
force acting on the tip (offsets added for clarity) exhibits a 
distinguished stick-slip movement at the step edge for down- 
and upward scans. Due to the relative increase of the bar- 
rier height and slope with load, the lateral force needed to 
pull the tip over the step edge increases for upward scans, 
(c) The tip-sample interaction potential V(xt) for the same 
three loads as in (b) , illustrating the increase of barrier height 
and steepness with load, (d) Plot of the frictional increase 
as calculated in the Prandtl-Tomlinson model for the poten- 
tial shown in (a). Parameters: Eo — 1.0 eV, rn = 0.45 nm, 
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1 kg, c x = 5.0 N/m and j x = -y z = 2 v 'c x m x . 



the force of the spring is large enough to pull it over 
the barrier. The force needed for this process, however, 
will typically be larger than on the terraces. After this 
first jump, the tip is temporarily stuck in the minimum 
at the bottom of the step edge. A second larger-than- 
usual jump is needed to disengage the tip, which occurs 
again when the condition Eq. © is fulfilled. Finally, the 
tip jumps out of the shallow minimum nested within the 
rising edge of the potential. 

For an upward scan [right to left in Fig.|4jc)], the tip 
first makes a "regular" jump from the last minimum of 
the periodic potential of the terrace into a shallow min- 
imum located within the potential's falling edge (first 
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solid arrow). From there, it takes a smaller-than-usual 
jump into the well developed minimum at the bottom of 
the step edge (second solid arrow). Pulling the tip now 
from this minimum over the barrier (third solid arrow) 
requires the application of a very large lateral force com- 
pared to typical forces experienced on the terraces. As a 
consequence, the tip subsequently jumps over a distance 
spanning several atomic unit cells, as can be seen from 
Fig.Hb). 

From the analysis, it is evident that the frictional forces 
are markedly different for downwards and upwards scans. 
Because the tip has a much higher barrier as well as a 
much steeper rise (slope of the potential) to overcome 
during the upward scan compared to the downward scan, 
lateral forces are naturally much higher for scanning up- 
wards. Also, we find that the relative height between 
the last minimum before the barrier and the maximum 
if approached from the side of the upper terrace does 
not significantly change with load, leading to a nearly 
load-independent frictional increase at the step edge. In 
contrast, both the relative barrier height as well as the 
slope increase if the barrier is approached from the side 
of the lower terrace, which leads to a linear dependence 
of the frictional increase on load. Fig.rjjd) illustrates this 
behavior, where the frictional increase has been plotted 
for upwards and downward scans. Note that these re- 
sults agree qualitatively very well with the experimental 
data of Figs. [I] and [2l suggesting that the experimentally 
observed load dependence is likely to be caused by the 
effects described in this paper. In this context, we would 
like to mention that similar results are obtained for a 
wide range of parameters as well as by representing the 
tip-sample interaction potential as a sum over Morse po- 
tentials. 

In summary, we presented experiments and simulations 
analyzing the load dependence of atomic-scale friction 
at surface steps. Experimentally, a direction-dependence 
has been found, where the contribution of the frictional 
forces due to the presence of the step edge increases 
linearly with load for upwards scans while it is load- 
independent for downward scans. By introducing a mod- 
ified Prandtl-Tomlinson model that includes an explicit 
description of the tip-sample interaction at surface steps, 
a theoretical basis for this behavior has been found. 
Finally, generalization of the above principles to other 
types of surface defects (vacancies, grain boundaries, 
etc.) where the atomic coordination will be temporarily 
altered suggests that this effect might very well dominate 
the macroscopic friction experienced on many materials. 

It is a pleasure to thank Andre Schirmeisen, Lars 
Jansen, Harald Fuchs (University of Minister) and Volker 
Saile (Forschungszentrum Karlsruhe) for discussions and 
their continuous support. This work was financially sup- 
ported by the BMBF (Grant No. 03N8704), the Na- 
tional Science Foundation (Grant No. MRSEC DMR 
0520495) , and the Petroleum Research Fund of the Amer- 



ican Chemical Society (Grant No. PRF 42259-AC5). 



Electronic address: hendrik.hoelscher@imt.fzk.de 
] D. Dowson, History of Tribology (Longman Inc., New 
York, 1979). 

] C. M. Mate, G. M. McClelland, R. Erlandsson, and 

S. Chiang, Phys. Rev. Lett. 59, 1942 (1987). 
■] B. Bhushan, J. N. Israelachvili, and U. Landmann, Na- 
ture 374, 607 (1995). 
J R. W. Carpick and M. Salmeron, Chem. Rev. 97, 1163 
(1997). 

] M. Urbakh, J. Klafter, D. Gourdon, and J. Israelachvili, 

Nature 430, 525 (2004). 
■] H. Holscher, A. Schirmeisen, and U. D. Schwarz, Phil. 

Trans. R. Soc. A (2008). 
] R. Liithi at al, J. Vac. Sci. Technol. B 14, 1280 (1996). 
■] S. Morita, S. Fujisawa, and Y. Sugawara, Surf. Sci. Rep. 
23, 1 (1996). 

'] J. Kerssemakers and J. T. M. De Hosson, Appl. Phys. 

Lett. 67, 347 (1995). 
i] H. Holscher, U. D. Schwarz, O. Zworner, and R. Wiesen- 

danger, Phys. Rev. B 57, 2477 (1998). 
] H. Holscher, et al, Phys. Rev. B 59, 1661 (1999). 
] A. Socoliuc, R. Bennewitz, E. Gnecco, and E. Meyer, 

Phys. Rev. Lett. 92, 134301 (2004). 
,] S. N. Medyanik, W. K. Liu, I.-H. Sung, and R. W. 

Carpick, Phys. Rev. Lett. 97, 136106 (2006). 
:] L. Prandtl, Z. Angew. Math. Mech. 8, 85 (1928). 
] G. A. Tomlinson, Philos. Mag. S. 7 7, 905 (1929). 
i] G. M. McClelland, in Adhesion and Friction, edited by 
M. Grunze and H. J. Kreuzer (Springer, Heidelberg, Ger- 
many, 1989), pp. 1-16. 
'] J. S. Helman, W. Baltensperger, and J. A. Holyst, Phys. 

Rev. B 49, 3831 (1994). 
.] E. Gnecco et al, Phys. Rev. Lett. 84, 1172 (2000). 
'] E. Riedo, E. Gnecco, R. Bennewitz, E. Meyer, and H. 

Brune, Phys. Rev. Lett. 91, 084502 (2003). 
i] A. Schirmeisen, L. Jansen, and H. Fuchs, Phys. Rev. B 

71, 245403 (2005). 
] G. Meyer and N. M. Amer, Appl. Phys. Lett. 57, 2089 
(1990). 

] J.-A. Ruan and B. Bhushan, J. Appl. Phys. 76, 8117 
(1994). 

■] E. Weilandt, A. Menck, and O. Marti, Surf. Interface 

Anal. 23, 428 (1995). 
:] E. Meyer et al, J. Vac. Sci. Technol. B 14, 1285 (1996). 
.] T. Midler et al, Phys. Rev. Lett. 79, 5066 (1997). 
■] We use the term "atomic-scale surface step" to refer to 
steps of one atom, two atoms, or another "low" number 
of atoms in height. Such steps are abundant on most 
crystalline surfaces. For an analysis of the lateral force 
signal observed at "large" surface steps (defined as steps 
where the tip radius becomes comparable or smaller than 
the height of the step) see, e.g., S. Sundararajan and B. 
Bhushan, J. Appl. Phys. 88, 4825 (2000), or I.-H. Sung, 
H.-S. Lee, and D.-E. Kim, Wear 254, 1019 (2003). 
[27] G. Ehrlich, J. Chem. Phys. 44, 1050 (1966). 
[28] R. L. Schwoebel and E. J. Shipsey, J. Appl. Phys. 37, 
3682 (1966). 

[29] U. D. Schwarz, P. Koster, and R. Wiesendanger, Rev. 



[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 
[8] 

[9] 
[10] 

[11] 
[12] 

[13 

[14] 
[15] 
[16] 

[17; 

[18] 
[19] 

[20] 

[21] 

[22] 

[23] 

[24] 
[25] 
[26] 



■5 



Sci. Instrum. 67, 2560 (1996). Sci. 375, 395 (1997). 

[30] H. Holscher, U. D. Schwarz, and R. Wiesendanger, Surf. 



